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The observed galaxy power spectrum acquires relativistic corrections from lightcone effects, and 
these corrections grow on very large scales. Future galaxy surveys in optical, infrared and radio 
bands will probe increasingly large wavelength modes and reach higher redshifts. In order to exploit 
the new data on large scales, an accurate analysis requires inclusion of the relativistic effects. This 
is especially the case for primordial non-Gaussianity and for extending tests of dark energy models 
to horizon scales. Here we investigate the latter, focusing on models where the dark energy interacts 
non-gravitationally with dark matter. Interaction in the dark sector can also lead to large-scale 
deviations in the power spectrum. If the relativistic effects are ignored, the imprint of interacting 
dark energy will be incorrectly identified and thus lead to a bias in constraints on interacting dark 
energy on very large scales. 

PACS numbers: 


I. INTRODUCTION 


There is no observational evidence yet for the exis¬ 
tence of non-gravitational interaction between dark en¬ 
ergy (DE) and dark matter (DM), nor is there guidance 
from fundamental theory as to the possible forms of such 
an interaction. However, there is also no consistent ex¬ 
planation for standard (non-interacting) DE from funda¬ 
mental theory. And interacting dark energy (IDE) is not 
ruled out by current observations. (For a range of work 
on IDE, see (T) 80] .) 

Various observations have been used to place con¬ 
straints on IDE models. In nearly all cases, the IDE 
models are compatible with current observations, pro¬ 
vided that the interaction is not too strong. A few of 
the models result in non-adiabatic large-scale instabilities 
(see, e.g. [40j [42]), and hence may be ruled out. How¬ 
ever, several of these models can be ‘fixed’ by adjusting 
the DE equation of state (see, e.g. nmmunu). 

As far as we are aware, there has not up to now been an 
analysis of the IDE imprints on structure formation that 
takes into account the general relativistic (GR) effects 
on the galaxy power spectrum. These GR effects arise 
from observing galaxies on the past lightcone, and they 
become significant on very large scales [HTj— |10lj . Fu¬ 
ture galaxy surveys will cover wide sky areas and reach 
high redshifts, and thus begin to probe scales beyond 
the the equality scale, and approaching the Hubble scale 
(at higher redshift z). In order to analyse observations 
on these scales, we need to use the correct theoretical 
model - i.e ., including the GR effects which deviate from 
the Newtonian approximation that is accurate on smaller 
scales. 

Future high-volume galaxy surveys will allow us to: 
(i) extend tests of DE and modified gravity models - and 
general relativity itself - to horizon scales; (ii) measure 
the primordial non-Gaussian signal in the galaxy power 
spectrum at higher precision levels than the CMB. GR 


effects have been discussed in the case of (i) by [951 M [98] 
and in the case of (ii) by [86j [881EU I1Q2H1Q4] . Recently, 
|105j investigated the degeneracy between the large-scale 
imprint of primordial non-Gaussianity and the imprint 
of some IDE models, but without considering the GR 
corrections to the power spectrum. 

Here we extend the work on DE models by investi¬ 
gating the GR effects on the galaxy power spectrum in 
IDE, assuming primordial Gaussianity. We normalise 
each IDE model to give the same matter density param¬ 
eter O m o and Hubble constant Hq as its corresponding 
standard non-interacting model. Then scale-dependent 
deviations from standard behaviour are clearly isolated 
on large scales at z = 0. For z > 0 however, the power 
spectra do not match on small scales. We show these 
behaviours in subsection IIV Bl 

The rest of the paper is organised as follows. We de¬ 
scribe the background IDE models in Section [TT] and the 
perturbed IDE models in Section III We investigate and 
discuss the very large scale behaviour of the power spec¬ 
trum, including GR corrections, in Section IV Then in 
Section [V] we conclude. 


II. BACKGROUND UNIVERSE WITH IDE 

In the background universe dominated by (cold) DM 
(.A = m) and DE (A = x), the interaction is defined via 
energy density transfer rates Qa , given by 

p' A + 3PL(1 + wa)pa = (1) 

where an overbar denotes background, a prime denotes 
derivative with respect to conformal time 77, and TL = 
a'/a is the conformal Hubble rate. The equation of state 
parameters are w rn = pm/pm = 0 and w x = p x /p x . 
The conservation of the total energy-momentum tensor, 
TT = 0? then implies 

Qm — Qx • 


( 2 ) 
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The case Q x > 0 corresponds to energy density transfer 
from DM to DE, and vice versa for Q x < 0. We may 
define effective equation of state parameters 

_ a QA /q\ 

WA,eff = WA~ • (3) 

3 np A 

where 

Pa + 3H(1 + WA,efi)PA = o, (4) 

so that w A ,eff encode the deviations from the standard 
evolution of the dark sector energy densities. 

The background field equations are determined by the 
total energy-momentum tensor X^, and thus do not 

explicitly contain the interaction: 


For non-interacting DE, the relativistic Poisson equation 
([8| then becomes V 2 4> = |H 2 J2a 

Note that the right hand side of the Poisson equa¬ 
tion © remains the same for both interacting and 
non-interacting DE, so that is always determined by 
^2 a ^a[^a — 3H(1 + w A )V A \• In the case of interact¬ 
ing DE, we have from ( [lT] ) that S A — 3H(1 + w A )V A = 
Au - aQ A V A /p A , thus 


V 2 $ 




Aa 


CiQ A 
PA 



( 12 ) 


A. Energy-momentum transfer 4-vectors 


■W 2 = ^(A, 


■Px), H' = --(l+3wA)n 2 , (5) 


where Q A = Pa/P and p is the total background energy 
density. (Note that p A implicitly contain the interac¬ 
tion.) We neglect baryons for simplicity - this does not 
affect the subject of our investigation. 


We discuss our various choices of Q A in Section III 


Generally, an interacting system is defined by the 
energy-momentum balance equations, given by nano 

no 


VArr = Q^ E^ = 0 ’ ( 13 ) 

A 

where Q A are the energy-momentum transfer 4-vectors 


III. PERTURBED UNIVERSE WITH IDE 


Qa ~ Qau» + F Al Qa — Qa + $Qa, u hF a — o. (14) 


The perturbed Friedmann-Robertson-Walker (FRW) 
metric in the conformal Newtonian gauge, with vanishing 
anisotropic stress, is given by 

ds 2 = a 2 [-(1 + 2 §)dri 2 + (1 - 24>)dx 2 ] , (6) 

where is the gravitational potential. The perturbed 
field equations do not contain the interaction explicitly 
when written in terms of Newtonian-gauge quantities, 
since they are determined by the total perturbed energy- 
momentum tensor a dT A u : 

& + m = ~n 2 ^A(i + w A )v A , (7) 

A 

vH = ln 2 Y J ^A[SA-2,n{i + w A )v A }. (8) 

A 

Here S A = dp a/Pa and the velocity potentials V A give 
the dark sector peculiar velocities, where the 4-velocities 
are 

<4 = a -1 (1 - &V A ), = a -1 (1 - &V) . (9) 

The total 4-velocity is u M and V is the total velocity 
potential, given by 

^1 + ^A^A^j V = ^ a (1 + w a)Va- ( 10 ) 

A A 

It is convenient to use the comoving overdensities 
A A = 5 A + ^V A = 5 A -3H(l + w A , eS )V A . (11) 

PA 


The energy density transfer rate Q A and the momentum 
density transfer rate F A are both relative to u 11 . In first- 
order perturbations, we have 

F a = a~ l (0, #f A ) , (15) 

where f A is the momentum density transfer potential. 
Then ([9|, ( fl()| ) and ( fl4| ) imply that 

Qo = — cl [Qa(1 + 4>) + 6Q a \ , 

Qt = adi [QaV + f A ] . (16) 

B. General perturbed balance equations 

By considering all species as perfect fluids, the per¬ 
turbed energy-momentum tensor of species A is given by 

STF = (6 p A + Sp A )u^ A u u A 

+ (pA + Pa) [Su a Ua + u A 5u A \ 

+ dp A r u +P A dg^, (17) 

where Su A and 5g are the perturbations in the 4- 
velocity and the metric tensor, respectively. The pressure 
perturbation Sp A is given by m 

5p A = c 2 A 6p A + 2,n{c 2 aA - c 2 a ) (1 + W A ,efl) p A V A , (18) 

where c 2 sA = (5p A /5p A )r es ttranie is the physical sound 
speed squared, and c 2 A = Pa/Pa * s th e square of the 
adiabatic sound speed. 
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Given (13)— (p~5]) , ( |17| ) and ([18]) we obtain the energy 
and the momentum 


balance equations: 


- 37-LwaAa + (1 + wa) V 2 Va 

- \u 2 {i + w a ) Y, n B (i + w B )\v A - v B ] 

B^A 

= 4[< 5Q a -MJa\ 

Pa 

+ n^fA A+Q ' {V - VA)] 

+ S jf{M(VA-V) + QAv A 

- 'i+-AMaA (19) 

L 1 + w A -I > 

V' A + HV a = -$ - 

1 + w A 

aQ A ( V - V A ) + af A , onA 

+ (1 + w a )pa ‘ ( j 

Here c srn = 0 and we take c sa , = 1 (which is the value for 
quintessence). 


C. Particular IDE models 

We model DE as a fluid with constant equation of state 
w x and consider two interactions. We use reCDM to de¬ 
note the non-interacting case and iw CDM the interacting 
case. 

In these interaction models, we assume that the trans¬ 
fer 4-vectors Q ^ run parallel to the DE 4-velocity: 

Qx =Qx<= -Qm- (21) 

This means that there is zero momentum transfer in the 
DE rest frame, which is the case for example in the mod¬ 
els of [9j [13j [T5j [25j [7TJ [76] [80]. From ( [21] ) , it follows 
that the momentum density transfer rates (relative to 
u become 


f x = Qx(v x - v) = -f m . 


( 22 ) 


For transfer 4-vectors of the form (|2l|, the balance equa¬ 
tions (19) - (20) lead to 


a ' m - ln 2 n x (i + w x )[v m -v x ] + v 2 v m 


3 U(V m - v x ) 


i 


aQ 

m 1 aQ rri 

31~ipm _ Pm 

T ~ \$Qrn T QmYm Qm^rn] 
Pm 


V' 


nvm = -$ - 


dQrr t 

Pm 


[V m ~ V x \ , 


(23) 

(24) 


for DM, and 


a ; - 3Hw x a x --H 2 n m (i + wx)iv x -v m 
= -(1 + w x )V 2 V x + 4 [SQx + Q' x Vx] 

Px 

J2 


1 + 


V' + nv x 


1 +w x 
-<I> 


Px 


l + W x 


(25) 

(26) 


for DE. (We have left c sx unspecified for generality, but 
we set it to 1 for numerical solutions.) 

To fully specify an IDE model, we need to define the 
Qa , which we choose as follows. 


Model 1: Q x oc p x 


We use a transfer rate [40 1 149 ] [105 ] . 


Qx — r p x — Tp x (l + 5 X ) — —Qm > (27) 

where T is a universal constant (i.e. it is fixed under 
perturbations). In the case T < 0, this corresponds to 
decay of DE into DM. 

Then from Q and p7| ), it follows that 

Ql = aYp x [-(1 + S x + <f>), diV x } = -Q™. (28) 


By (22), 


fx = rp x (V x -V) = -f m . (29) 


Model 2: Q x oc ©p x 


It is common in the literature to use energy den¬ 
sity transfer rates of the form Q x = £,a~ 1 Hp Xl where 
£ = const, i.e. to use a transfer rate proportional to the 
Hubble rate, rather than a constant rate T as Model 1 
(27). The motivation for this choice is that the back¬ 


ground energy conservation equations are easily solved. 
However, the problem is that for the perturbed model, 
the Hubble rate T-L is typically not perturbed. 

We resolve this problem by using instead the self- 
consistent transfer rate 


Qx — 3 £ p *0 j 


0 = VX- 


(30) 


In the background, this reduces to the usual form, but in 
the perturbed universe we pick up additional perturba¬ 
tions of the expansion rate: 


© = 3a -1 n - ($' + H$) + \v 2 V 

o 


(31) 
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FIG. 1: Evolution of the IDE effective equation of state parameters for the w CDM equation of state parameters 

w x — —0.8 (left panel) and w x = —1.1 (right panel). Solid lines correspond to Model 1 (27) while dashed lines correspond to 


Model 2 (|30j) , and the w x line denotes T = 0 = £. 


This leads to 


Q x = a~ 1 Oip x [l + S x - 4> - — (3$' - V 2 V) J = -Q m . 

(32) 

Then (jo|, ( [21] ) , (31) and (32) imply that 


Ql = m x 


-i ~ 5x + m ~ ’ diVx 


• (33) 


By (22) 


fx = a~^np x (V x -V) = -f n 


(34) 


For both models, the range of w x is restricted by sta¬ 
bility requirements 


w x > —1 for £, T > 0; w x < —1 for £, T <0. (35) 


These two cases correspond to different energy transfer 
directions, by (27) and (30): 

DM DE for f, T > 0; DE —» DM for f, T < 0. (36) 


D. Background evolution of the models 


We evolved the background equations from around de¬ 
coupling, ad = 10 —3 , until today ao = 1. The background 
initial conditions were chosen so that the matter density 
parameter o and the Hubble constant Hq are the same 
as in the non-interacting models. We used f2 m o = 0.24, 
H 0 = 73 kms _1 Mpc _1 . 

It is easy to understand the behaviour of the effective 
equation of state parameter w XjG s in Fig. ITJ For Model 1 

(07 \ o 11 ... — _o "p / ( A ~ ^ i □ o 


(27), w X:e ff = w x — aF/(3'H), where al~L~ is a positive 


growing quantity. Hence w XjG s gradually decreases for 
T > 0 (w x = —0.8) and increases for F < 0 (w x = —1.1), 
with |T| determining the strength of the interaction (solid 
lines). However, for Model 2 (30) we have w XjG ff = w x — 
£/3, which is a constant for constant w x , shown in the 
dashed lines. 

Figure [2] shows the evolution of the matter density 
parameters and Hubble rates, compared to the non¬ 
interacting case. When \Q X \ is small (|£|, \F / Hq\ < 0.1), 
w Xie ff is only slightly less (greater) than w x = —0.8 
(w x = —1.1). This implies weaker (stronger) DE effects: 
the background matter density for iwCDM becomes en¬ 
hanced (suppressed) relative to reCDM for w x = —0.8 
(w x = —1.1). However, when the transfer rate is higher 
m,\r/H 0 \ >0.4), w x , e ff is much smaller (bigger) than 
w x = —0.8 (w x = —1.1) and hence iwCDM has sur¬ 
plus (less) matter relative to reCDM for w x = —0.8 
(w x = -1.1). 

Notice the distinct separation between the ratios of the 
Hubble rates. To understand this, we know that during 
matter domination, when « 1 (Q x <C 1), the ratio 
is constant. In this work we fixed the background ini¬ 
tial conditions in reCDM and let those in iw CDM vary 
with each value of F or £ so as to recover the same val¬ 
ues of U m0 and H 0 as in reCDM. As T or £ vary, the 
initial conditions in iw CDM change, enough to amount 
to significantly differing initial amplitudes of the Hubble 
rates. The ratio does not evolve until DE domination, 
converging at a = 1 by our normalization. 


IV. THE LARGE-SCALE POWER SPECTRUM 

We probe the large-scale structure of the late-time Uni¬ 
verse by relating the perturbations at late epochs to the 
primordial potential via growth functions and a transfer 
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FIG. 2: Ratios of the matter density parameters (left) and Hubble rates (right) for zieCDM relative to those of rcCDM: with 
w x = —0.8 (top panels) and w x — —1.1 (bottom panels). Line styles are as in Fig. [l] 


function (see subsection IVA), from which the matter 
power spectrum is computed. 


A. Linear growth functions 


Here is the gravitational potential at photon-baryon 
decoupling. It is related to the gravitational potential 
growth function by [98] : 


$d(k) 


D^a) D*{k,a d )=a d , 

^ p (k)T(k), 

, e m0 /_fc_\ (n_4)/2 

D* 0 (k) \H 0 ) 


(37) 

(38) 

(39) 


where D$> 0 = D$(k, 1), is the primordial poten¬ 
tial and T(k) is the transfer function which accounts 
for perturbation evolution through radiation domination 
until radiation-matter transition. The constant A = 


h\/2 tt5h /(3i?g /2 ) is the primordial amplitude of curva¬ 
ture perturbations, we adopt 5h = 5.6xl0 -5 (see [53] ) 
which is the scalar amplitude at horizon crossing and 
n = 0.96 is the scalar spectral index. 

The growth function D m of the comoving matter over¬ 
density A m describes the growth of (linear) matter per¬ 
turbations after radiation-matter equality via [98 

A m (k,a) =--f—^D m (k,a)^ d (k). (40) 

o\LmO 

The matter velocity growth function Dy m is defined by 
V m (k,a) = - 2 tj 2 Dv m (fc, a)^rf(fc), (41) 

Oi L rn QllQ 

where we have used 0 and assumed matter domination 
and Einstein de Sitter regime (i.e. with T' = 0). Note 
that Djn implicitly contains the interaction: 

D m (k,a) = D° m (k,a) + ^AAD Vm (k,a), (42) 
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FIG. 3: The ratio of the matter overdensity and gravitational 
potential growth functions, at ao = 1 (or z = 0), with w x = 
—0.8. Solid lines correspond to Model 1 (27) and dashed lines 
to Model 2 (30). The ACDM case (dashed black line) and the 
Hubble horizon (solid black line) are also shown. 


where D ^ is the growth function in the standard non¬ 
interacting scenario. Equation (42) is obtained by using 
(37) and (40) in the Poisson equation (12) - and assumed 
matter domination. 


Note that D<$> and Dv m in (37) and (41) retain the same 
definitions in both IDE and non-IDE scenarios - i.e. they 
do not explicitly contain the density transfer rate Q m - In 
standard ACDM, the growth functions become 


D m (k,a) = D<s>(k,a) = D(a), Dy m (fc,a) = 


dD(a) 

dr] 


.(43) 


However, for dynamical DE (w x ^ —1), D m = D<$> only 
holds true on sub-Hubble scales k T-L. 

The DE velocity growth function is related to that of 
matter by 


D Vx (k, a) 


(^0 


[1 - fi m (a)] [1 + w x (a)\ dr) 
-D Vm {k,a) 


d 

— D<s>(k, n) 


(44) 


where we have used (37) and (41) in (|7|) and assuming 
Einstein de Sitter regime - i.e. with 4>' = 0. 

The DE overdensity growth function is 


Dx (k. cl) — 


ifl) 

1 (n) - 

^Pm (&) 


D$(k, a) Djyiiyk^ci) 

Dv x (k, a) — D Vm {k, a) j ,(45) 


where we used ( |37| ), © and ( |40| ) in the Poisson equa¬ 
tion (12), and the fact that Q x Jp x = — D m (2 m /(O x p m ). 
In ACDM, by (43), D x = 0 = Dy x . 


B. Relativistic effects in the galaxy power 
spectrum 


Using cosmological perturbation theory, we can de¬ 
scribe the evolution and distribution of matter density 
perturbation in the Universe. However, in reality the 
matter perturbation distribution is not directly observ¬ 
able - only objects such as galaxies, whose distribution 
traces that of the matter, are observable. The large-scale 
and scale-independent galaxy bias b is usually defined by 
S g = bS m . On sub-Hubble scales, different gauge choices 
for S m agree, but on near- and super-Hubble scales, they 
disagree. The scale-independent bias needs to be defined 
physically, in the rest-frame of DM and galaxies (which 
coincide on large scales). This leads to the following def¬ 
inition, valid on all linear scales (and assuming Gaussian 
primordial perturbations) [85j 8j)J |88| : 

A g (fc, a) = b(a) A m (fc, a). (46) 

Furthermore, we do not observe in real (x) space, but 
in redshift space, leading to the Kaiser redshift-space dis¬ 
tortion term: 


Ag bs (fc, /i, a) = [b(a) + f(k,a)/j, 2 ]A m (k,a), 
where 

f = Dvm 
J UD m 


(47) 

(48) 


which reduces to the growth rate of matter overdensity 
in non-interacting DE models, p = —n • k jk and n is the 
unit spatial vector in the direction of the photon geodesic 
x^(\) from source to observer. Here A is the affine param¬ 
eter along the photon geodesic, increasing from source to 
observer. (Note that / in (48) is different from f m and 
fx in (|22|.) 

The Kaiser term is a flat-sky and sub-Hubble approx¬ 
imation to the full redshift-space distortion, which in¬ 
cludes further velocity and Sachs-Wolfe type terms. In 
addition to the redshift distortion, there are other rela¬ 
tivistic effects from observing galaxies on the past light- 
cone. These include a contribution from weak lens- 
ing convergence, which can be significant on sub-Hubble 
scales at higher redshifts. In addition there are addi¬ 
tional Sachs-Wolfe and Doppler terms, and integrated 
Sachs-Wolfe (ISW) and time-delay terms [8TJ |82[ [8H l85l 
18811931 [106]. 


If we want an accurate analysis that includes near- and 
super-Hubble scales, we should use the galaxy overden¬ 
sity that is observed on the lightcone, including all GR ef¬ 
fects. This observed overdensity is automatically gauge- 
invariant. Here we will neglect the integrated terms and 
use a flat-sky approximation, generalizing the ACDM 
form given in [88] : 


A f s {k, ii, a) = A f(k, )i, a) + A ° R (fc, n, a), (49) 


where A( ul is the standard term - given by ( |47| ). Note 
that in ([49]) we do not have a priori the time-delay, ISW 
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FIG. 4: Observed galaxy power spectrum Pg bs (solid lines) and the standard power spectrum PJ td (dashed lines) along the line 
of sight (= 1), at z = 0: for w x = —0.8 {top); and for w x = —1.1 (bottom). The corresponding ratios of the power spectra are 
given on the right panels. 


and weak lensing integrated terms, and 


A gr = 


A 


k 2 


■ i/iB 


H 


(50) 


This form arises by using the field and conservation equa¬ 
tions to relate velocity and potential to overdensity. (See 
e.g. [93, 97, 98 for the full GR expression including in¬ 
tegrated terms.) The coefficients in (50) are given in the 
ACDM case by [88]. We generalize their expressions for 
the case of IDE: 


A = (3 b e ) f 


B 


+ 2 


-2 


2 n 2 D m , 

(1 -Q) a 2 a /D*\ 

rU + D$da\ a Jl 

(1 - Q) , aQ 
rU 


4Q - b e - 1 + —, 


W_ 

n 2 


ki. (si) 

a 


+ ^.( i _A i )] / 

rlprn ^ -Dvl ' - 




(52) 


Here r is the comoving radial distance at the observed 
galaxy, Q is the magnification bias, and b e is the galaxy 
‘evolution bias’, giving the evolution of source counts (see 

[HI M [88]). 


The interaction enters B through the last term in (52). 
This term arises from the perturbed Euler equation (24), 
which comes in via 


^A.(n-V m ) = ^n-V l m + ^n%(n-V m ), (53) 

where the total derivative is taken along the photon 
geodesic, in the direction from source to observer. The 
Q m term is absent in ACDM [8TJ [82j [84| [S3 Ell and in 
non-interacting DE models [98]. It would also be absent 
in IDE models with Q^ parallel to the matter 4-velocity 
since in these models the DM follows geodesics and 
the perturbed Euler equation is the same as for non¬ 
interacting DE. 


From (49)—(52), we obtain the power spectrum P° bs of 
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FIG. 5: Ratios of the observed galaxy power spectrum P° bs to the standard power spectrum P| td along the line of sight (/jl = 1), 
at z = 1: for w x — —0.8 (left) and w x — —1.1 (right). 


the observed galaxy over density in the conformal Newto¬ 
nian gauge (we give only the real part): 


:>obs 


n 2 


(b + ffJ?) +2(6 + /M 2 )^-p- 


nj 4 nj 2 

+ A 2 — - u 2 B 2 — 
^ ^ k 4 ^ ° k 2 ’ 


(54) 


The standard power spectrum ( i.e . the Kaiser approxi¬ 
mation) is given by 


P| td (fc, fi, a) = [6(a) + f(k, a)/j 2 ] 2 P m (k, a). (55) 


The matter power spectrum is obtained from <|40|) and 
(|37b-(|39l» as [98j 


Pm {k^ ck) 


9 A 2 

507r 3 Py 


k n T(k) 2 


D m (k,a) I 2 
D* 0 (k) 


(56) 


We note that the terms in the last line of ( [54] ) correspond 
to the auto-correlations of Aj^ R , and the last term in 
the first line corresponds to the cross-correlation between 
A| td and A™. 

To compute the growth functions, we employ adia¬ 
batic initial conditions (see Appendix 0 . For the goal 
of this work, it is reasonable to assume a constant co¬ 
moving galaxy number density so that b e = 0. We also 
set (henceforth) the galaxy bias to b = 1 for simplicity. 
We considered the case with the magnification bias 


2 = 1 , 


(57) 


which corresponds to intensity mapping of the neutral hy¬ 
drogen (HI) 21cm emission line [95l 1981 [TOT] . In this case, 
the lensing convergence and time delay terms drop out of 
Ag bs , and the remaining integrated term, the ISW term, 
typically makes a negligible contribution. This gives a 


strong justification for our neglect of the integrated terms 
in @-([52]). 

Equation (56) shows that the induced changes in P m 
arising from the interaction will be imprinted via the ra¬ 
tio D m /D<s>o. This ratio is exactly unity on all scales in 
ACDM while for dynamical DE it tends to 1 (by nor¬ 
malization) only on small scales k P, where the DE 
perturbations are negligible. This is illustrated in Fig. [3] 
using iwCDM with w x = —0.8 as an example. In this 
case, DM loses energy to DE, causing suppression of the 
matter growth function, which shows up in the matter 
power spectrum. Note that in the absence of interaction 
(as in the standard cosmologies), the clustering of DE 
causes large-scale suppression in the matter power m. 
Moreover, the growth functions show that Model 1 (27) 
is relatively more sensitive to the values of the interaction 
parameter, compared to Model 2 (30). 

For w x = —1.1 in iw CDM, the corresponding plots for 
Dm /simply have all the curves reflected across that 
of ACDM. 

In Fig. [4] we present the line-of-sight GR corrected 
galaxy power spectrum Pg bs and the standard power 
spectrum P| td in the Kaiser approximation, for the 
iwCDM models at z = 0. On sub-Hubble scales, (54) 
shows that 


pobs 

g 


(1 + /) 2 for k^>%, 


(58) 


i.e. P° bs 


pstd 

g • 


We also see that ireCDM gives a large-scale boost in 
galaxy power for both T, £ > 0 (w x > — 1) and T, £ < 0 
(w x < — 1). As remarked above, there is less sensitiv¬ 
ity to the values of |£|, so that Model 2 ( |3Q| ) predicts 
relatively larger amplitudes on super-horizon scales than 
does Model 1 (27). 

The large-scale boost in galaxy power in Fig. [4] for 
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k (Mpc 1 ) 



k (Mpc x ) 



— r/H 0 =o=(; 

— r/H 0 =-o.i 

— r/H 0 =-0.4 
— £=— 0.1 

— 4 = -0.4 


-3 



k (Mpc 1 ) 


FIG. 6: Ratios of P° bs , which is P° bs with the correlation between A| td and A^ R subtracted, to the standard power spectrum 
P| td , along the line of sight (fi = 1). Top panels have z = 0, for w x = —0.8 (top left) and w x = —1.1 (top right). Similarly for 
z = 1 in the bottom panels. 


rcCDM arises purely from GR effects. The smaller boost 
(relative suppression) in power for iw CDM with r, £ > o 
(top left) comes from the fact that DM loses energy to 
DE, so that the higher the DM rate of loss of energy 
(i.e. larger |T|, |£|), the more the power suppression. For 
T, £ < 0 (bottom left), where DE loses energy to DM, 
larger values of |T| and |£| give more boost relative to 
wCDM. 

Here, GR corrections in the galaxy power spectrum re¬ 
sult in large-scale galaxy power enhancement. This im¬ 
plies that if we ignore the GR effects, i.e. if we do not 
subtract them in order to isolate the IDE effects, then 
we arrive at an incorrect estimate of the imprint of IDE 
on very large scales. 

In Fig. [4] (right panels), we also show the correspond¬ 
ing ratios of the total galaxy power spectrum relative to 
the Kaiser approximation, at z m 0. The ratios show 
that, for T, £ > 0 (top right), GR effects are suppressed 


relative to reCDM by the interaction. This is consis¬ 
tent with previous explanations above. Conversely, for 
T, £ < 0 (bottom right), GR effects are enhanced rela¬ 
tive to reCDM. 


However, note that the total galaxy power spectrum 
contains not only the individual contributions of the stan¬ 
dard Kaiser redshift-space distortion and GR terms, but 
also their cross-correlation. This cross-correlation makes 
a positive contribution at low z, and a negative contri¬ 
bution at high z - for the given magnification bias (57). 
Hence, the ratios shown (at z = 0) contain positive con¬ 
tributions of this cross correlation term, as well as the 
auto-correlation of the GR corrections. The ratios at 
z = 0 are thus enhanced on horizon scales. 


The GR effects in our IDE models, where we use 
Gaussian primordial perturbations, show a similar be¬ 
haviour to the effects of primordial non-Gaussianity in 
non-interacting DE models (with /nl > 0). The degener- 
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acy between GR effects and primordial non-Gaussianity 
in the ACDM model has been investigated by [86] [551 
l94l H02H104j . It was recently shown by [105] that IDE 
effects (in the case where GR effects are neglected) can 
be degenerate with primordial non-Gaussianity. 

In Fig. [5j we show the ratios of the observed galaxy 
power spectrum to the Kaiser approximation at higher 
redshift, z = 1. For higher values of the interaction pa¬ 
rameters, the case with T, £ > 0 (left) shows lower large- 
scale GR effects in comparison to the case with T, £ < 0 
(right). For the magnification bias ( [57] ) (corresponding to 
HI intensity mapping), the observed line-of-sight power 
spectrum falls to zero for both non-interacting and inter¬ 
acting DE. 

At higher z the IDE effects are weaker, since the ef¬ 
fects of DE in general are weaker at earlier times. By 
contrast, the GR effects are typically stronger at higher 
z - but with A^ R having negative amplitude. Hence its 
correlation with the Kaiser term gives negative contribu¬ 
tion in the power spectrum, thereby gradually reducing 
galaxy power on horizon scales. For completeness, we 
illustrate this phenomenon in Fig.[6j which shows the ra¬ 
tio of P° bs , which is Pg bs with the correlation between 
A| td and A^ R subtracted, to P| td . In the top panels, 
z = 0, and z = 1 in the bottom panels. On the left, 
w x > — 1 and T, £ > 0, and on the right, w x < —1 and 
T, £ < 0. By comparing the top left and top right panels, 
with the top right and the bottom right panels of Fig. [4j 
respectively, we see that the correlation between the GR 
and the Kaiser terms is positive at low z. Similarly, by 
comparing the bottom left and the bottom right panels, 
with the left and right panels of Fig. [5j respectively, we 
see that the cross correlation is negative and of larger 
amplitude at z = 1. 


V. CONCLUSION 


Future wide and deep-field surveys may be able to dis¬ 
entangle any possible IDE effects from GR effects by com¬ 
paring the observed power at low and high z. Detecting 
super-Hubble effects will be challenging because of cos¬ 
mic variance. However, if the multi-tracer method [108] 
can be applied, cosmic variance can be reduced enough 
for detection of these effects m- 
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Appendix A: Adiabatic initial conditions 

We use the Einstein de Sitter initial condition 4>'(a^) = 
0, given that Q x (a,d) <C 1. Adiabatic initial conditions are 
imposed by the vanishing of the relative entropy pertur¬ 
bation S xm , 


Sxmip'd) — 0 ? S X m — 


1 T ^Uc.eff 1 T 


(Al) 


and by the equality of velocities, 


V x (ad ) = Vm(ad ). (A2) 

Then (1 T = (1 T VJx,efi (^d)) A m (tt^). 

Together with ( [I2| and ([7]), this leads to the initial DM 
and DE fluctuations given by 


We investigated GR effects in the observed galaxy 
power spectrum in two IDE models, comparing with 
the corresponding standard non-interacting DE scenar¬ 
ios. We focused on the case of magnification bias Q = 1, 
corresponding to HI intensity mapping, and normalized 
the IDE power spectra to those of their non-interacting 
DE counterparts on small scales at today, i.e. by requir¬ 
ing that they have the same U m o and Hq. This isolates 
the deviations arising from GR effects and IDE on very 
large scales. 

We find that if the GR effects are disregarded, i.e. if 
they are not subtracted in order to isolate the IDE effects, 
then we arrive at an incorrect estimate of the imprint 
of IDE on horizon scales. This could lead to a bias in 
constraints on IDE on the given scales. 

We also found that at low z, the correlation between 
the GR term and the (standard) Kaiser redshift-space 
distortion term has a positive contribution in the galaxy 
power spectrum, while at high z, this term gives a nega¬ 
tive contribution that grows with increasing z. 


A m (k) = 
A x (k) = 


2,h (1 ^7n,eff) 


(i “i - 

_ ~ 2 fc 2 (1 + Wx,eff) _ 

3(1 “1“ ff 

v x {k) = V m (k) = ~ 2 - -Mk), 

on. (1 + \l x w x ) 


MV, (A3) 
*<*(*), (A4) 
(A5) 


where w me ff and w xx ,u are given by ([ 3 ]) and <!>,/(A') is 
given by (38). 
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